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Strongly correlated many-body systems show various exciting phe- 
nomena in condensed matter physics such as high-temperature su- 
perconductivity and colossal magnetoresistance. Recently, strongly 
correlated phases could also be studied in ultracold quantum gases 
possessing analogies to solid-state physics, but moreover exhibit- 
ing new systems such as Fermi-Bose mixtures and magnetic quan- 
tum phases with high spin values. Particularly interesting systems 
here are quantum gases in optical lattices with fully tunable lattice 
and atomic interaction parameters. While in this context several 
concepts and ideas have already been studied theoretically and 
experimentally, there is still great demand for new detection tech- 
niques to explore these complex phases in detail. 

Here we report on measurements of a fully momentum-resolved 
excitation spectrum of a quantum gas in an optical lattice by 
means of Bragg spectroscopy. The bandstructure is measured with 
high resolution at several lattice depths. Interaction effects are 
identified and systematically studied varying density and excita- 
tion fraction. 

The understanding of excitations is crucial for any physical system. Es- 
pecially in many-body physics it is important to characterize a system's 
response in different regimes to obtain information about the underlying 
phenomena. Techniques such as angle-resolved photoemission spectroscopy 
(ARPES) measuring the spectral function, or inelastic neutron scattering 
giving the dynamical structure factor are examples which have greatly con- 
tributed to advances in solid-state physics. Measuring these quantities in 
quantum gases in optical lattices is a very promising step towards further 
connecting these model systems with their condensed matter counterparts. 

Important progress has been made in this direction in lattice-free geome- 
tries, e.g. the spectral function has recently been measured in trapped Fermi 



gases[l] and the dynamical structure factor 5(A:,li;) has been successfully 
studied by optical Bragg spectroscopy in free and trapped Bose-Einstein 
condensates [2, 3, 4, 5, 6] as well as trapped Fermi gases[7]. 

Measurements of excitations in optical lattices have been performed by e.g. 
a potential gradient[8] and parametric excitations[9, 10] without momen- 
tum resolution. Our experimental setup allows for the first time a freely 
variable momentum transfer through the first Brillouin zone and therefore 
gives access to the whole excitation spectrum. This provides great oppor- 
tunities studying the properties of phases realizable in these systems such 
as a Mott-insulator, a Tonks- Girardeau gas or supersolid phases as recently 
proposed[ll, 12, 13, 14, 15, 16, 18, 17]. Up to date, measurements in optical 
lattices employing Bragg spectroscopy have been restricted to a single[19] 
or two different momenta[20, 21]. 

Momentum-resolved Bragg spectroscopy 
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Figure 1: Schematic representation of the Bragg process in real and momentum space. 

a Bragg process in real space: Two laser beams with momentum hki and energy licoi crossing 
at an angle O excite an ensemble of atoms caught in an optical lattice, b Bragg process in 
momentum space: The figure shows the excitation spectrum of a square lattice. The Bragg pro- 
cess transfers energy and momentum via a two-photon process and populates a different state 
thereby giving information on the excitation spectrum. Here an excitation in [l,l]-direction 
close to the Brillouin zone edge is depicted. 

Bragg spectroscopy in quantum gases is based on a two-photon process 
which directly transfers energy and momentum to an ensemble of atoms. As 
the experimental situation is fundamentally different to similar techniques 
in condensed-matter physics, a few aspects are elaborated for clarity. Two 
phase-coherent laser beams with wavevectors ki and k2 and frequencies uji 



2 



and L02 intersect on the sample as depicted in Figure la. Due to energy and 
momentum conservation, this process resonantly couples initial and final 
states with momentum difference ?i/cBragg and energy difference h5 = h{L02 — 
uji). Figure lb visualizes the Bragg process in momentum space. In contrast 
to e.g. ARPES, the momentum is directly imprinted on the sample and is 
freely variable by changing the angle G between the two beams: 

hkBragg = Kh - ~ 2hksm (0/2) • (Cfcj - 4 J (1) 

with e/c^ , the unit vectors in the propagation direction of beam 1 and 2 
respectively, and ~ |/i;2| = /c for small 6. As the energy transfer h5 can be 
tuned independently by changing the difference of the two laser frequencies, 
we can actively scan the excitation spectrum across the Brillouin zone with 
full access to all bands. 

In order to study the excitation spectrum of an ultracold bosonic cloud of 
atoms in the optical lattice, we perform Bragg spectroscopy for more 
than 20 different Bragg laser angles @ and several lattice depths Vq. For 
each /sBragg and each lattice depth, we vary the energy difference 6 between 
the beams and measure resonance spectra as depicted in Figure 3a. 

Bandstructure in optical lattices 

Exposing atoms to a standing light wave creates a periodic potential for neu- 
tral atoms. Due to the periodicity, a bandstructure evolves and the particles 

(n) 

can be described by Bloch functions (pj, with quasimomentum k, band index 
n and energy in the case of non-interacting atoms. Quantum gases 

usually show short-range interactions characterized by a s-wave scattering 
length oo leading to an interacting many-body system. This is described in 
good approximation by the Bose-Hubbard model 

H = - J ^ a)^aRt + ^ X] "-^(^^ ~ ^) + X] ^R^R 
{R,R') R R 

with J the hopping matrix element and U the on-site interaction[22] . The 
operator fiR counts the number of atoms at lattice site R, er is the external 
potential and {R, R!) restricts the tunneling a^a^/ to nearest neighbours. 

In order to study such a system, we measure its excitation spectrum for 
different parameters and evaluate the dependencies on crucial values. As 
a central result. Figure 2 shows for the first time the detailed excitation 
spectrum of a quantum gas in an optical lattice ranging from small k up to 
the edge of the Brillouin zone for different lattice depths Vq. For our mea- 
surements we have chosen the momentum transfer along the [l,l]-direction 
in a square lattice, also called the nodal direction from F to M. Figure 2a 
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Figure 2: Excitation spectra, a Excitation energies with respect to momentum transfer fcsragg 
for harmonic trapping potential(A) and lattice depths Vo = 3Er( ), 7Er(B) and llEr(^): 
The solid lines show calculations including interactions. Statistical errors always lie within the 
data points, b Excitation energies versus momentum transfer at higher lattice depths Vo = 7Er 
and 11 Er! The dashed lines show the single-particle bandstructure while the solid lines show 
the results of a Bogoliubov-approximated Bose-Hubbard calculation (see text). Error bars give 
the statistical uncertainty in energy, the uncertainty in momentum transfer lies within the data 
points. 



displays the resonance positions of the lowest excitations corresponding to 
the first band for several lattice depths. The band is well resolved and its 
decreasing width with increasing lattice depth is visible. A spectrum of a 
harmonically trapped quantum gas is shown for comparison. Note that com- 
pared to electrons, the excitations of bosonic atoms in optical lattices show 
a distinct phonon-like behaviour in the low-momentum regime, embodied by 
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Figure 3: Resonance spectra and corresponding bandstructure. a Relative population in 
the central peak versus laser detuning 5 at ^Bragg ~ for lattice depths Vo — 3Er(#), 

7Er(B): The two main peaks of each spectrum correspond to the first and second band, 
respectively. Fit functions are given by a sine function as the system's response to a square 
pulse, b Resonance position (first and second band) with respect to momentum transfer feg^ggg 
along the [l,l]-direction: Dotted lines correspond to a single-particle bandstructure calculation. 
Solid lines show the excitation energies considering interactions based on a numerical solution 
of the Bogoliubov-de Gennes equation for the second band and following equation (3) for the 
first band. Dashed lines show a numerical solution of non-interacting particles with an attached 
external confinement (see text). Error bars give the statistical error in energy, the uncertainty 
in momentum transfer always lies within the data points. 



a linear increase in the excitation energy, as has been studied in free Bose- 
Einstein condensates [4]. In this regime, the lattice reduces the sound velocity 
with higher lattice depth as expected[23]. Close to the Brillouin zone edge. 
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the influence of the lattice fundamentally changes the excitation spectrum 
resulting in the opening of a bandgap which grows with increasing lattice 
depth. Figure 3b shows this exemplarily for the lattice depths Vq = SE^ 
and 7E]-, with Er = ti^k\/ {2m) the recoil energy and the wavevector of 
the lattice laser. The influence of the lattice at the Brillouin zone edge also 
reflects significantly the dynamics of the system. In the absorption images 
which map the population in momentum space two peaks appear at inter- 
mediate momentum transfers. They correspond to the initial state at /c = 
and the final state at A; = A^Bragg- Getting close to the Brillouin zone edge, a 
third peak at /c = /cBragg — G, with G a reciprocal lattice vector, evolves and 
becomes more and more pronounced. We attribute this to Bragg refiection 
at the zone edge, i.e. Bloch oscillations [24]. 

Figure 2b shows a comparison of the experimental data to both the theory of 
a non-interacting single-particle and a mean-field approximation of a Bose- 
Hubbard model for the lattice depths Vq = 7Ei. and 11 Ef. As the dashed 
lines show the results of a single-particle bandstructure calculation, it is 
obvious that interaction strongly modifies the excitation spectrum even in 
the superfiuid regime. It is now especially interesting to check the validity 
of a mean-field approach for quantum gases in optical lattices within a wide 
range of parameters. To start, we compare our data with a Bogoliubov-type 
approximation in the Bose-Hubbard framework following Burnett et al.[25] 
which gives an analytical expression for the excitation energies 



hwj^ = W4Jsin^ 



ka 

Y 



o I ka 

2nU + 4Jsm^ ( y 



(3) 



with a the lattice spacing. Since the momentum transfer is in [l,l]-direction 
with kx = ky, we obtain excitation frequencies cok^^ky = + ^ky- These 
are depicted as solid lines in Figure 2b showing the expected resonance 
frequencies for our experimental parameters at the corresponding lattice 
depths. While the interaction effects are still relatively small for Vq = 3Er, 
they increase substantially with larger lattice depth owing to the rise in the 
term 2nU in equation (3). The deeper the lattice, the stronger the atoms are 
confined locally at the lattice sites resulting in a larger on-site interaction U. 
In addition, a deeper lattice also leads to a stronger harmonic confinement 
and thus to an increase in the average atom number per site n. The good 
agreement with the experimental data shows the applicability of a mean-field 
treatment in this regime. 

An inherent advantage of optical Bragg spectroscopy is the easy access to 
higher bands. In Figure 3b the bandstructure showing both the first and 
the second band at lattice depths Vb = 3Er and 7Er is displayed. The 
dotted lines show the results of a single-particle bandstructure calculation 
whereas the solid line is obtained by a numerical solution of the Bogoliubov- 
de Gennes equation which includes interaction on a mean-field level[26]. 
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Interaction effects in tlie second band are ratlier small due to the small 
overlap of initial and final state wavefunctions. This is consistent with a 
harmonic oscillator approximation of a single lattice site and the symmetry 
of the corresponding wavefunctions of the ground and first excited state. 
However, especially for a shallow lattice of Vq = 3Er, both calculations 
cannot reproduce the data as we observe much lower excitation energies at 
small momentum transfer. We attribute this shift to the influence of the 
confinement, since, considering zero-point and excitation energy, the second 
band lies energetically above the lattice potential. We incorporate this by 
solving a single-particle problem in a finite lattice with an attached harmonic 
confinement. The solution is depicted by the dashed line in Figure 3b and 
shows good agreement with the experimental data. In the case of the deeper 
lattice of Vq = 7Er, the second band lies almost completely within the lattice 
potential and the effect is drastically reduced. 
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Figure 4: Shift of resonance position with density. Resonance position versus transversal 
lattice depth: The lattice depth in 2D is always Vb = 5 Er, and the measurements are performed 
at three momenta fegragg = 0.64^ (■), ^Bragg = 0-71^( ) and fegragg = 0.93f;(A). The 
inset show/s resonance frequency versus atom number for a 2D(A) and a 3D lattice(#) w/ith 
transversal lattice depth Vr = 5Er at S = 45° and otherwise identical parameters. Error bars 
show the statistical uncertainty in energy. 



Influence of interactions 

As we have discussed above, the first band shows large interaction effects. 
In order to quantitatively characterize the influencing factors of interaction, 
we systematically study the shift of the resonance position with respect to 
variation of density, on-site interaction JJ as well as excitation fraction in 
the Bragg process. 
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We investigate the influence of the density n and on-site interaction U on 
the excitation energy by adding a third lattice axis, which is orthogonal to 
the Bragg process. In contrast to the tight binding electron orbitals in solid- 
state physics, the wavefunctions of a quantum gases in optical lattices are 
expected to separate in spatial dimensions and are, in a first approximation, 
composed of harmonic oscillator orbitals. First we study this assumption 
by varying the atom number in a 2D and 3D lattice at otherwise identi- 
cal trapping and Bragg parameters and measure the shift of the excitation 
energy. The inset of Figure 4 shows the corresponding data points and fits 
based on equation (3). While we consider J as fixed for a given lattice depth 
in the lattice plane, the term 2nU is described within the Thomas-Fermi 
limit [27]. The rise in excitation energy due to higher particle number and 
due to the transversal lattice axes is very well reproduced within our model. 
Thus we can estimate the relative increase in nU from a 2D lattice to a 3D 
optical lattice. Even more important, it proofs that within the experimen- 
tal accuracy the transversal lattice depth does not infiuence the tunneling 
J in the 2D lattice plane which is consistent with a separation ansatz for 
the wavefunctions. This allows us to independently change the parameter 
nU via the transversal lattice axis. Figure 4 shows the resonance frequency 
versus the transversal lattice depth for three different momenta, while the 
lattice depth in the 2D square lattice is fixed at Vq = 5Er. The increase in 
excitation energy with rising transversal confinement is clearly visible. Based 
on equation (3), for a certain lattice depth, we only change the momentum 
k in this type of experiment. This allows the extraction of the evolution of 
2nU and an independent measure of the experimentally given confinement 
with good accuracy. 

We conclude that the transversal lattice axis does not influence the excita- 
tion spectrum apart from changing density and on-site interaction. Further- 
more, we can extract relevant parameters and resolve shifts in the excitation 
energy on the order of a few hundred Hz due to either density or on-site inter- 
action which have to be taken into account for spectroscopic measurements 
as well as state preparation [28]. 

By varying the pulse area, we can control the fraction of transferred atoms. 
For the linear response regime small excited state fractions are required, 
while larger transfer ratios are important in state preparation and studies 
of the final state. In combination with the total number of atoms, we can 
change the density in the initial state ipi and also in the final state ijjf. This 
provides further insight into the nature of the final state and gives access to 
the differential interaction between ■i/'j and ipf. Figure 5 shows the particle 
number dependence of the resonance position for different pulse areas in the 
magnetic trap. We vary the excited state fraction between 3% and nearly 
100%. The excitation energy is lowest in the case of a full population transfer 
realized by a vr-pulse, where we also see the gradual increase with particle 
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Figure 5: Resonance position versus atom number with respect to the excitation fraction. 

Excitation energy versus atom number at ^Bragg ~ 0.71^ for harmonic trapping potential: 
Measurements are performed at excitation fractions of 75%(A), 45%(B), 20%( ), 6%(4) 
and 3%(A), respectively. The inset shows the corresponding data in a 2D lattice geometry at 
a lattice depth Vo = 7Er. Measurements are performed at excitation fractions of 80%(A), 
45%(A), 23%( ) and 8%(*), respectively. Error bars show the statistical uncertainty in en- 
ergy. 



number, i.e. density. Reduction of the excited state fraction leads to higher 
excitation energies. The effect saturates at excitation fractions of 10% to 
20% which we identify as the boundary of the hnear response regime for 
the free case. The inset of Figure 5 shows the corresponding measurements 
in the optical lattice exemplarily at a lattice depth Vq = TEf. We observe 
a similar behaviour, however, the size of the effect seems to be strongly 
reduced. Within our accuracy, saturation already sets in below a fraction 
of 25% given by the green circles. Controlling the fraction of atoms in the 
final state reveals the change in interaction energy within the initial state 
Uii, within the final state Ujj and between initial and final state Uif. Looking 
at the slope of the increase in excitation energy with respect to the atom 
number, we derive in this simple picture that interactions Uif and Uff are 
smaller than Uu considering the same number of particles. In the lattice, the 
difference between Ua and Uff decreases significantly, although corrections 
are still clearly visible. Slightly adjusting this technique by performing a 
state-changing Raman transfer and using a Feshbach resonance to tune Uif 
to zero would complement this information to determine the change of the 
wavefunction with momentum and offers further insight into the complex 
many-body state [29]. 

Determining the pulse area necessary for a full transfer in each spectroscopic 
measurement gives an estimate for the efficiency of the process and therefore 
the matrix element connecting the two states. The efficiency in the first band 




9 



increases with higher momentum transfer up to the edge of the first Brillouin 
zone as expected[ll]. Looking at Figure 3a ihustrates that first and second 
bands exhibit different peak heights in the spectrum at otherwise identical 
parameters. The transfer to the second band at this momentum transfer 
seems to be more efficient than that to the first band which differs from 
the above mentioned calculations by Menotti et al.[ll] where a reduced 
excitation efficiency is predicted. This aspect also shows the suitability of the 
method to determine the static structure factor which is a quantity that can 
be used to detect and distinguish specific phases such as a supersolid[17, 18]. 

In conclusion we have measured a fully momentum-resolved bandstructure 
of a bosonic quantum gas in an optical lattice. We could quantify many-body 
interaction effects with high resolution and compare to different theoretical 
models. At the same time, Bragg spectroscopy turns out to be a powerful 
tool to prepare states with a chosen momentum or even a superposition of 
several momenta. This work paves the way for further studies of quantum 
phases in optical lattices and is especially interesting for fermionic systems 
as well. 



Methods 

We routinely produce a EEC of ^'^Rb in a decompressed magnetic trap 
with a trap frequency of a; = 27r • 14 Hz and an average particle number of 
N ~ 150000. After ramping up the optical lattice in 2 or 3 dimensions in 
150 ms, we start the Bragg process after 10 ms of hold time. As shown in 
Figure la and lb, two beams detuned to an excited state by a frequency 
A ~ 5 GHz with a freely variable angle are used to realize a momentum 
transfer in [l,l]-direction. The momentum transfer is given in units of ^), 
with a the lattice spacing a = 515 nm. Thereby the value 1.0 corresponds to 
the Brillouin zone edge in the nodal direction. 

We employ two different experimental sequences in the low- and high-momentum 
regime which are described below. Detection of the atomic clouds is per- 
formed via absorption imaging after a time-of-flight (TOF) of ixoF = 20 ms. 

We used two different methods to determine the resonance position: coher- 
ent transfer and energy transfer. In the case of coherent transfer, the Bragg 
pulse is applied for a time = 1 ms and lattice as well as magnetic trap are 
switched off immediately afterwards. This results in projecting the quasimo- 
mentum states on real momentum states: In the case of no Bragg transfer, 
these are the main cloud centered around zero momentum and lattice peaks 
having momentum of an inverse lattice vector G around 2hkL^x,y Applying 
the Bragg pulse transfers momentum hkBragg to the excited atoms. During 
TOF the atoms move according to their momentum and separate clouds 
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become visible after a certain time. An absorption image is taken and the 
fraction of transferred atoms at feeragg as well as the reduction of atoms at 
A; = is determined versus the frequency difference of the beams resulting 
in a resonance spectrum (see Figure 3a) . 

Due to the small momentum transfer at lower angles, the clouds do not 
separate any more on experimental timescales and we determine the reso- 
nance condition by energy transfer to the system. The Bragg pulse excites 
atoms and imparts a certain amount of energy into the system. We allow 
the system to rethermalize before switching off lattice as well as magnetic 
trap and measure the width of the atomic cloud or the condensate fraction 
after TOF. 

We checked for low, intermediate and high momenta and for different lattice 
depths that both methods, coherent transfer and energy transfer, yield the 
same results. 

Acknowledgements 

We gratefully acknowledge valuable discussions with C. Becker, U. Biss- 
bort, W. Hofstetter, K. Rachor and K. Bongs. We also thank the Deutsche 
Forschungsgemeinschaft DFG for funding within the Forschergruppe FOR801 

Competing financial interest 

The authors declare no competing financial interests. 

References 

[1] Stewart, J.T., Gaebler, J. P. & Jin, D.S. Using photoemission spec- 
troscopy to probe a strongly interacting Fermi gas. Nature, 454, 744- 
747 (2008). 

[2] Kozuma, M. et al. Coherent Splitting of Bose-Einstein Condensed 
Atoms with Optically Induced Bragg Diffraction. Phys. Rev. Lett., 
82, 871-875 (1999). 

[3] Stenger, J. et al. Bragg Spectroscopy of a Bose-Einstein Condensate. 
Phys. Rev. Lett, 82, 4569-4573 (1999). 

[4] Stamper-Kurn, D.M. et al. Excitations of Phonons in a Bose-Einstein 
Condensate by Light Scattering. Phys. Rev. Lett, 83, 2876-2879 
(1999). 



11 



[5] Ozeri, R., Katz, N., Steinhauer, J. &; Davidson, N. Colloquium: Bulk 
Bogoliubov excitations in a Bose- Einstein condensate. Rev. Mod. Phys., 
77, 187-205 (2005). 

[6] Papp, S.B. et al. Bragg Spectroscopy of a Strongly Interacting ^^Rb 
Bose-Einstein Condensate. Phys. Rev. Lett, 101, 135301 (2008). 

[7] Veeravalli, C, Kuhnle, E., Dyke, P. & Vale, C.J. Bragg Spectroscopy of 
a Strongly Interacting Fermi Gas Phys. Rev. Lett, 101, 250403 (2008). 

[8] Greiner, M., Mandel, O., Esslinger, T., Hansch, T.W. & Bloch, I. Quan- 
tum Phase Transition from a Superfluid to a Mott Insulator in a Gas 
of Ultracold Atoms. Nature, 415, 39-44 (2002). 

[9] Stoferle, T., Moritz, H., Schorl, C, Kohl, M. & Esslinger, T. Transition 
from a Strongly Interacting ID Superfluid to a Mott Insulator. Phys. 
Rev. Lett, 92, 130403 (2004). 

[10] Schorl, C, Stoferle, T., Moritz, H., Kohl, M. & Esslinger, T. Excitations 
of a Superfluid in a Three-Dimensional Optical Lattice. Phys. Rev. 
Lett, 93, 240402 (2004). 

[11] Menotti, C, Kramer, M., Pitaevskii, L. &; Stringari, S. Dynamic struc- 
ture factor of a Bose-Einstein condensate in a one-dimensional optical 
lattice. Phys. Rev. A, 67, 053609 (2003). 

[12] Roth, R. & Burnett, K. Dynamic structure factor of ultracold Bose and 
Fermi gases in optical lattices. J. Phys. B: At. Mol. Opt. Phys., 37, 
3893-3907 (2004). 

[13] Hofstetter, W., Cirac, J.I., Zoller, P., Demler, E. & Lukin, M.D. High- 
Temperature Superfluidity of Fermionic Atoms in Optical Lattices. 
Phys. Rev. Lett, 89, 220407 (2002). 

[14] van Oosten, D., Dickerscheid, D.B.M., Farid, B., van der Straten, P. &; 
Stoof, H.T.C. Inelastic light scattering from a Mott insulator. Phys. 
Rev. A, 71, 021601 (2005). 

[15] Rey, A.M., Blakie, P.B., Pupillo, G., WilHams, C.J. & Clark, C.W. 
Bragg spectroscopy of ultracold atoms loaded in an optical lattice. Phys. 
Rev. A, 72, 023407 (2005). 

[16] Pupillo, G., Rey, A.M. & Batrouni, G.G. Bragg spectroscopy of trapped 
one-dimensional strongly interacting bosons in optical lattices: Probing 
the cake structure. Phys. Rev. A, 74, 013601 (2006). 

[17] Mathey, L., Danshita, I. & Clark, C.W. Creating a supersolid in one- 
dimensional Bose mixtures. Phys. Rev. A, 79, 011602 (2009). 



12 



[18] Titvinidze, I., Snoek, M. &; Hofstetter, W. Supersolid Bose- Fermi Mix- 
tures in Optical Lattices. Phys. Rev. Lett, 100, 100401 (2008). 



[19] Du, X., Wan, S., Yesilada, E., Ryu, C. & Heinzen, D.J. Bragg spec- 
troscopy of a superfluid Bose-Hubbard gas. http://arxiv.org/pdf/ 
0704.2623 (2007). 

[20] Fabbri, N. et al. Excitations of Bose-Einstein condensates in a one- 
dimensional periodic potential, http://arxiv.org/pdf/0901.1805 
(2009). 

[21] Clement, D. et al. Exploring correlated ID Bose gases from the 
superfluid to the Mott-insulator state by inelastic light scattering, 
http : //arXiv . org/0812 . 4530V1 (2008). 

[22] Jaksch, D., Bruder, C, Cirac, J.I., Gardiner, C.W. & ZoUer, P. Cold 
Bosonic Atoms in Optical Lattices. Phys. Rev. Lett., 81, 3108-3111 
(1998). 

[23] Kramer, M., Menotti, C, Pitaevskii, L. & Stringari, S. Bose-Einstein 
condensates in ID optical lattices: compressibility, Bloch bands and 
elementary excitations. Eur. Phys. J. D, 27, 247 (2003) 

[24] Dahan, M.B., Peik, E., Reichel, J., Castin, Y. & Salomon, C. Bloch 
Oscillations of Atoms in an Optical Potential. Phys. Rev. Lett., 76, 
4508 (1996) 

[25] Burnett, K., Edwards, M., Clark, C.W. & Shotter, M. The Bogohubov 
approach to number squeezing of atoms in an optical lattice. J. Phys. 
B: At. Mol. Opt. Phys., 35, 1671-1678 (2002). 

[26] Berg-S0rensen, K. &; M0lmer, K. Bose-Einstein condensates in spatially 
periodic potentials. Phys. Rev. A, 58, 1480-1484 (1998). 

[27] Pitaevskii, L. & Stringari, S. Bose-Einstein Condensation. (Oxford 
University Press, Oxford, 2003). 

[28] Miiller, T., FoUing, S., Widera, A. & Bloch, L State Preparation and 
Dynamics of Ultracold Atoms in Higher Lattice Orbitals. Phys. Rev. 
Lett, 99, 200405 (2007). 

[29] Dao, T.-L., Georges, A., Dalibard, J., Salomon, C. & Carusotto, I. Mea- 
suring the One-Particle Excitations of Ultracold Fermionic Atoms by 
Stimulated Raman Spectroscopy. Phys. Rev. Lett, 98, 240402 (2008). 



13 



